arXiv:astro-ph/9701133v2 25 Jan 1997 


THE STRUCTURE OF THE X-RAY EMITTING GAS IN THE 


HYDRA-A CLUSTER OF GALAXIES 

Y. Ikebe^ K. Makishima^, H. Ezawa^, Y. Fukazawa^, M. Hirayama^’^, H. Honda^, 
Y. Ishisaki^’^, K. Kikuchi"^, H. Kubo^’^, T. Murakami^, T. Ohashi^, T. Takahashi^, 

and K. Yamashita^ 

ikebe@postman.riken.go.jp: maxima@miranda.phys.s.n-tokyo.ac.jp: 
ezawa@tkyosfl.phys.s.n-tokyo.ac.jp: fnkazawa@miranda.phys.s.n-tokyo.ac.jp: 
hirayama@tkyosfl.phys.s.u-tokyo.ac.jp: honda@astro.isas.ac.jp: 

ishisaki@phys.metro-n.ac.jp: kiknchi@phys.metro-n.ac.jp: 
hknbo@miranda.phys.s.u-tokyo.ac.jp: murakami@astro.isas.ac.jp: 
ohashi@phys.metro-u.ac.jp: takahasi@tkynext.astro.isas.ac.jp: 
yamasit a@sat io. phys. nagoya- u. ac. j p 

Received_; accepted_ 


^Cosmic Radiation Laboratory, The Institute of Physical and Chemical Research 
(RIKEN), 2-1 Hirosawa, Wako-shi, Saitama 351-01, JAPAN 

^Department of Physics, University of Tokyo, 3-1, Hongo 7-chome, Bunkyo-ku, Tokyo 
113, JAPAN 

^The Institute of Space and Astronautical Science, 1-1, Yoshinodai 3-chome, Sagamihara- 
shi, Kanagawa 229, JAPAN 

“^Department of Physics, Tokyo Metropolitan University, 1-1, Minamioosawa, Hachioji- 
shi, Tokyo 192-03, JAPAN 

^Department of Physics, Nagoya University, Furoucho, Chikusa-ku, Nagoya-shi, Aichi 


464-01, JAPAN 





- 2 - 

ABSTRACT 

The temperature and abundance structure in the intracluster medium (ICM) 
of the Hydra-A cluster of galaxies is studied with ASCA and ROSAT. The 
effect of the large extended outskirts in the point-spread function of the X-Ray 
Telescope on ASCA is included in this analysis. In the X-ray brightness prohle, 
the strong central excess above a single /5-model, identihed in the Einstein and 
ROSAT data, is also found in the harder energy band (>4 keV). A simultaneous 
£t of five annular spectra taken with the GIS instrument shows a radial 
distribution of the temperature and metal abundance. A significant central 
enhancement in the abundance distribution is found, while the temperature 
profile suggests that the ICM is approximately isothermal with the temperature 
of ~3.5 keV. The ROiSCTPSPC spectrum in the central 1'.5 region indicates a 
significantly lower temperature than the GIS result. A joint analysis of the GIS 
and PSPC data reveals that the spectra can be described by a two temperature 
model as well as by a cooling flow model. In both cases, the hot phase gas with 
the temperature of ~3.5 keV occupies more than 90% of the total emission 
measure within 1'.5 from the cluster center. The estimated mass of the cooler 
(0.5-0.7 keV) component is ~2-6xl0® Mq, which is comparable to the mass 
of hot halos seen in non-cD ellipticals. The cooling flow model gives the mass 
deposition rate of 60 ± 30 Mq yr“^, an order of magnitude lower than the 
previous estimation. 


Subject headings: galaxies:clusters:individual(Hydra-A) — X-rays:galaxies 
X-rays: spectra 
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1. INTRODUCTION 


X-ray imaging studies with the Einstein observatory and the ROSAT satellite revealed 
that many clusters exhibit central concentrations in their X-ray brightness prohles. The 
central excess emission is often interpreted as a result of a cooling flow, a thermal instability 
occurring at the densest part of cluster (see |habian 1994] for a review). Edge et ah (1992) 
reported that about two thirds of all known clusters have such structure, suggesting that 
this structure is a fairly common characteristic among clusters. Spectroscopic studies with 
Einstein and ROSAT found that there are cooler gas components with temperature n-u 10^ K 
in cluster centers (e.g. [Uanizares et al. 1979| ). However, the previous spectroscopy is limited 
to the soft X-ray band (~0.1-4 keV). ASCA has provided the hrst opportunity to resolve 
X-ray spectra spatially in the wider 0.5-10 keV energy range. 

In this paper we report the ASCA observation of the Hydra-A cluster of galaxies 
{z = 0.0522). Its X-ray luminosity measured with the Einstein Observatory was 4.1x10^"^ 


ergs s " (0.5-4.5 keV) ([David et al. 19901) , and a consistent value of 4.8 x 10^“^ ergs s ^ in 
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the same energy band was indicated from the Ginga observation which measured the X-ray 
spectrum in the 2-10 keV band ( j'hsuru 1991| ). These luminosities are among the largest of 
poor clusters ( jKriss, Cioffi, fc Canizares 1983] ; [Tsuru 19'M|) . The average X-ray temperature 
of ~4 keV, as measured with Ginga, is also rather high for a poor cluster. 

The central region of this cluster is of particular interest. The cD galaxy 3C218 = 
Hydra-A is a strong radio source, with a very complex radio morphology and the highest 
Faraday rotation ever measured from a radio galaxy (pjlkers fc Simkin 1983| ; [Kato et "aTj 
|1987| ; [laylor et al. 199'0|) . Furthermore, the X-ray surface brightness prohle obtained with 
the Einstein IPC exhibits a large central excess deviating from an isothermal /d-model, 
suggesting the presence of a cooling flow with an estimated mass deposition rate of 600 ± 
120 Mq yr“^ ([David et al. 19901) . This is one of the largest values thus far attributed to the 
X-ray emitting gas in clusters of galaxies ([Edge et al. 1992[) . 

When we analyze the ASCA data, the complex response of the X-Ray Telescope (XRT; 
Serlemitsos et al. 1995[) onboard ASCA introduces considerable difficulty for spatially 
resolved spectroscopy. To confront with these difficulties, some analysis methods have been 
proposed ([Ikebe 1995[ ; [Markevitch et al. 1996[; [Churazov et al. 1996[) . The method employed 
in this paper is based on the similar idea in Markevitch et al (1996) but is performed 
with different implementation. Since the position resolution of ASCA is relatively poor, 
the sharp ROSAT image is also very helpful to analyze the ASCA data. Assisted by 
the ROSAT image, we measured the spatial distribution of the temperature and metal 
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abundance in the Hydra-A cluster from the ASCA observation. In section 2, we describe 
the ASCA observation and the data selection procedure. Section 3 gives the results using 
a conventional analysis method on the GIS and SIS data. Section 4 describes the analysis 
method developed here. In section 5, we analyze the X-ray brightness prohles in different 
energy bands. In section 6, we describe the temperature and abundance structure studied 
using the GIS spectra, and we also discuss the joint analysis of the ROSAT and ASCA 
data to study the temperature structure in the central region. We discuss the results and 
summarize them in section 7. Hq = 50 km s“^ Mpc“^ is assumed throughout. 


2. ASCA OBSERVATION AND DATA SELECTION 


ASCA ( [Tanaka, Inoue fc Holt 1994|) observed the Hydra-A cluster during the 
performance verihcation phase on 1993 May 28 and 29. The two GIS detectors (jOhashi et 
jal. 19961 ; [Makishima et ah 1996| ), GIS-S2 and GIS-S3, were operated in PH normal mode. 
The two SIS sensors, SIS-SO and SIS-Sl, were operated in 4GGD faint or bright mode, and 
central region of the Hydra-A cluster was detected in the middle of chip 1 and chip 2 of 
SIS-SO, and chip 3 and chip 0 of SIS-Sl. 


We have applied the standard data cleaning procedure (see pay et al. 19951) to the 
GIS and SIS data, and selected the data with the criteria that the cut off rigidity and 
elevation angle must be greater than 8 GeV/c and 5.0 degrees, respectively. The additional 
selection criterion that the elevation angle from the sunlit earth is greater than 20 degree 
was applied to the SIS data only. In the present analysis, we used 24.4 ksec of the SIS-SO 
data, 22.7 ksec of the SIS-Sl data and 26.3 ksec of the GIS data survived the selection. The 
X-ray intensity contour map from the summed GIS-S2 and S3 data is shown in Fig. |l|. As 
the background (cosmic X-ray background plus non X-ray background) of the GIS data, 
we utilized the data from the Lynx field, which was observed on 13-15 May 1993 when the 
spread discriminator in the onboard GPU of GIS was disabled as it was for the Hydra-A 
cluster (see [Makishima et al. 1996|) . For the SIS background, we used the data of several 
blank sky helds observed during the PV-phase. 
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3. ANNULAR SPECTRAL ANALYSIS 


For detailed study of the spatial variations in the spectra, energy spectra are often 
accumulated from different regions and are analyzed individually (see e.g. Phashi et ah 


19941 ; |Ikebe et ah 199'|) . However, this conventional analysis method is no longer usable 


for the ASCA data of clusters of galaxies. Although the ASCA XRT has a sharply peaked 
Point-Spread Function (PSF) with about 3 arcmin half-power diameter (HPD), the PSF 
has wide outskirts with signihcant flux extending over the whole GIS held of view of 52 
arcmin diameter (Serlemitsos et ah 1995|) . Therefore individual regions on the focal plane 
have inevitable contribution of the hux from other regions. Moreover, such a hux-mixing 
effect depends considerably on the X-ray energy because the PSF has wider outskirts at 
higher energies. In particular, this energy-dependent hux-mixing ehect is quite serious for 
the study of clusters of galaxies. Since their X-ray surface brightness decreases from center 
to periphery by several orders of magnitude, the outer regions of clusters inevitably have 
a considerable and energy-dependent hux from the inner, much brighter parts. Therefore 
the meaning of the individual spectra accumulated over specihed areas in the focal plane is 
very unclear. Figure demonstrates the mutual contribution of X-ray hux from annular 
sky regions to the annular regions on the focal plane in case of the Hydra-A cluster. 


Although not accurate, in this section we hrst perform model httings to individual 
annular spectra, as have been done in previous papers on the ASCA data. This is to 
illustrate the results when the energy-dependent hux-mixing ehect is not considered, and to 
compare them with the results from the new analysis method described in the subsequent 
sections. 


The GIS and SIS spectra have been accumulated in concentric annular regions centered 
on the cD galaxy, covering radius ranges of 0 — 1'.5 (arcmin), 1'.5 — 3', 3' — 5', 5' — 10', 
and 10' — 20'. The data of GIS-S2 and GIS-S3 were then combined after the correction 
for the gain diherence. For the SIS, data from four diherent chips in the same detector 
were combined after relative gain corrections. The SIS and GIS spectra thus obtained both 
exhibit strong emission lines from highly ionized ions. We htted annular spectra individually 
with a thermal bremsstrahlung model plus Gaussians modihed by photo-electric absorption, 
using the XSPEG package (ver. 9.0). The energy ranges used for the htting are 0.7-10 keV 
and 0.35-10 keV for the GIS and SIS, respectively. For all the spectral fittings, we used the 
XRT effective area which is calculated assuming that the point source is being observed with 
a circular detector of 6' radius size. Figure ^ shows the central 0 — 1'.5 spectra together with 
the best-fit model. At least three Gaussians at the rest-frame energies of 6.7 keV, ~2.0 keV 
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and ~1.0 keV are required, which may correspond to the He-like Fe-Ko,, He-and H-like 
Si-Kcj, and a mixture of He-like Fe-Lc^ and H-like Ne-K„ emission lines, respectively. There 
is also some indication of a line structure at ~1.5 keV, which may consist of He-like Fe-L^ 
and H-like Mg-Ko, emission lines. The continuum is well described by a single temperature 
model. When we performed the htting using only below 3 keV energy band or above 3 keV 
energy band, the best-£t temperatures were not different. The two temperature model did 
not improve the best-£t values at all. We illustrate the projected radial prohles of the 
temperature and the equivalent width (EW) of the Fe-K emission lines in Fig. Unless 
otherwise specihed, all errors are 90% confidence for one parameter of interest (Ay^ = 2.71). 
The two instruments have given consistent results, and the gradual increase towards the 
outside is seen in the temperature prohle. 


We also htted the spectra with the single-temperature thin-thermal plasma emission 
model by Raymond and Smith (1977; hereafter RS-model). The temperature, heavy element 
abundances, and normalization were varied as free parameters, while relative abundance 
ratios were hxed at the solar values given by Anders & Grevesse (1989). The projected 
radial prohles of temperature and abundance derived are shown in Fig. The temperature 
prohle obtained exhibits a gradual increase towards the outside and is consistent with the 
result from the hts with a thermal bremsstrahlung model plus Gaussians (Fig. §-(a)). The 
abundance values obtained with the SIS are signihcantly smaller than those obtained with 
the GIS (Fig. |-(b)). The same systematic diherence has been found in the data from 
various clusters observed with ASGA ( |Fukazawa 1996|) . This may be because the actual 
SIS energy resolution is slightly worse than the prediction of the response matrix, probably 
because of the chip summation. When we applied an alternative plasma emission model, 
Mewe-Kaastra model (MEKA-model; |Mewe, Gronenschild, fc van den Oord 1985t [Mewe 
fc van den Oord 1986| ; [Kaastra 19921) , we derived consistent results with those from the 
RS-model. Gonsistency between RS-model and MEKA-model results has been found from 
the ASGA data of other clusters with temperature of ~3-4 keV ([Mushotzky et al. 1996|) . 
As already pointed out, we cannot interpret the radial prohles presented in Fig. ^ and Fig. ^ 
as true projected structure on the sky. The gradual increase towards the outside seen in 
the temperature prohles is not a real temperature structure but is likely to be an artifact 
due to the XRT responses. Simulations assuming an isothermal distribution also produce a 
similar temperature prohle (I'Fakahashi et al. 1995|). 


4. ANALYSIS METHOD FOR THE EXTENDED SOURCES 
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4.1. Analysis Methods 


One way to analyze extended sources is to perform simultaneous fitting of the X-ray 
image and the spectra; that is, simultaneous htting of the spectra extracted from different 
regions, or similarly, simultaneous htting of the images sorted in different energy bands. We 
start with an X-ray source model which, in case of a cluster, describes spatial distributions 
of the X-ray emissivity, temperature, and abundance. A given model is converted to the 
simulated ASCA data through the convolution calculations with the instrument responses 
which include the XRT PSF and effective area, as well as the quantum efficiency and energy 
redistribution function of the GIS or SIS. The goodness of ht of the model is evaluated by 
comparing the simulated data with the observed data, and the results are iteratively fed 
back to the input model. 


The XRT-I-GIS PSFs used in our analysis were produced by interpolating the actually 
obtained X-ray images of Gyg X-1 with GIS. However, the Gyg X-1 is so bright that the SIS 
can not measure its brightness prohle precisely due to event pile-up; thus can not reproduce 
the XRT-I-SIS PSFs with the same method. Therefore, in this paper we analyzed only 
the GIS data with the analysis method described here. The detailed description of PSF 
calibration is in the Appendix. 


A technical problem in this analysis is how to reduce the calculation time. The image 
convolution, which is the most time consuming part, is mathematically expressed as an 
operation by the 2-dimensional matrix called the image response matrix with size as large 
as ~ 20000 X 20000 for the GIS full-resolution image. To reduce the calculation time, 
we have developed an analysis method using the Monte-Garlo technique ([lakahashi et al.| 
1995 ; Ikebe 1995f) which is often much faster than the huge matrix operation. We only 


need to generate events more numerous, say by a factor of 10, than the actual events. 
This Monte-Garlo method was employed for the analysis of the Fornax cluster ([Ikebe et al. 


|1996|) , the Gentaurus cluster ( jlkebe 19'95|) , 3A0336-I-098, A1795, MKW3s and PKS2354-35 
( Phashi et al. 1995] ; [Kikuchi et al. 19^) , Ophiucus cluster (jMatsuzawa et al. 1996|) , and 
Goma cluster ([Honda et al. 1996|) . 


4.2. Image Response Matrices 

An alternative method to reduce the calculation time is to reduce the size of the 
image response matrix. Due to the limited number of observed photons, we usually bin 































the image so that each image pixel has sufficient number of events that statistics can 
be employed. Moreover, an X-ray source model given is often so simple that the number 
of sky regions can be largely reduced. Thus, the dimension of the image response matrix 
can be reduced drastically. Since the Hydra-A cluster can be reasonably modeled with a 
circularly symmetrical model, as will be described in the following sections, we employed 
the matrix operation to generate simulated data rather than using the Monte-Carlo method. 
This method was hrst proposed by Markevitch et al (1996) and had been applied to the 
ASCA data analysis of A2163, A754, A2256, A2319, A665, Triangulum Australis, A3558, 
and AWM7 ( [Markevitch et al. 1996| ; [Henriksen fc Markevitch 1996| ; [Markevitch 1996] ; 
Markevitch, Sarazin, fc Irwin 199^ ; [Markevitch fc Vikhlinin 1997[) . Here we recall the 
formula of the matrix operation and dehne the function for the minimum htting. 

The model prediction is given by the form: 


MODEL%,= Y. 


sky Qsky 

energy energy ^ energy 


( 1 ) 


energy sky 

where IBM is the image response matrix, which is dehned as: 


TTDj^det, sky _ pOTpdet, sky pppsky 
liuvl energy JAOr energy ^ energy * 


( 2 ) 


Here, S is the initial cluster model on the sky which gives the X-ray surface brightness as a 
function of X-ray energy (expressed with the suffix energy) and sky position {sky ); EFF is 
the effective area including quantum efficiency of the focal plane instruments expressed as 
a function of X-ray energy and sky position; PSF is the XRT PSF which gives the photon 
distribution on the detector {det) depending on the X-ray energy and sky position; and 
BA'IF is the energy redistribution function which expresses pulse height {pi) distribution for 
an incident monochromatic-energy X-ray event. While both the ASCA data and the initial 
cluster model on the sky (S') have three dimensions, the htting is performed in the form of 
some energy-sorted brightness prohles, or in the form of some ring-sorted spectra. 


To evaluate the best-ht parameters of the assumed model, we dehned the function 


as: 


X' 


= EE 

det 


{DATA^' 


et 

pi 


B3D% 


MODEL'^'P] 


pi 




(3) 


where 

, + (“ X BGDXf + (6 X MODELS? ■ (4) 

Here DATA, BGD and MODEL are the values of the observed data counts, the normalized 
background, and the model, respectively, a^yj^j^det and are statistical errors 
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associated with respective quantities, while the latter two terms in eq. (4) represent 
systematic errors. The systematic error of the image response matrix is included in the 
value by adding 5% (6 = 0.05) of the model flux. We also introduced 10% of the background 
intensity estimated from the blank sky data in the systematic errors; that is a = 0.1. 

We have checked the consistency between the Monte-Carlo method and results from 
the image response matrices. For example, from a best £t model obtained by htting using 
the image response matrices, the Monte-Carlo simulation can reproduce the data which 
agree with the actual data with sufficient accuracy. 


5. RADIAL BRIGHTNESS PROFILE 

5.1. ASCA Data and Fitting Resnlts 

Using the analysis method described in § 4.2, we hrst analyzed the X-ray images taken 
with the CIS. The X-ray image has a good circular symmetry, and the azimuthally-averaged 
background-subtracted radial brightness prohles in different energy ranges are shown in 
Fig-B They are centered on the peak of the X-ray intensity and have 20 bins in total of 1 
arcmin width each. With ASCA, we were able to obtain the image in the 4-10 keV range 
for the hrst time, and these radial prohles are very similar among diherent energy bands. 

In order to perform the model htting to the radial prohle following the method 
described in § 4.2, we constructed image response matrices of a size 26 x 20 for each of 
the 201 energy bins. Each matrix represents hux contribution from 26 sky regions to the 
20 regions on the focal plane at a specihed energy. The 26 sky regions employed here are 8 
annular regions of O'.25 width for r < 2' and 18 annular regions of 1' width for r = 2' — 20', 
while the detector region was divided into 20 annular bins of 1' width each. Since the matrix 
used for the htting represents contribution only from sky regions within 20' radius, it is 
implicitly assumed that no emission is generated outside of the 20' radius. This assumption 
is consistent with the ROSAT image which will be described in the next section, and even 
when we assume the maximum radius to be larger than 20', the htting results do not change 
signihcantly. 

As the model brightness prohle, we employed an empirical (3 model expressed as: 



S(r) = So 


(5) 
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where r is the projected angular distance from the center, (3 represents the beta parameter 
and Tc is the core radius. We assume that the energy spectrum takes the same form in the 
entire cluster and is expressed with the RS-model modihed by the Galactic absorption. The 
temperature, heavy element abundance, and the hydrogen column density are assumed to 
be 3.4 keV, 0.5 solar and 6 x 10^° cm“^ respectively, as derived in § 3 from the spectral 
htting of the GIS data extracted from the inner most region. The free parameters are Tc, (3 
and Sq. 

With this single /3-model, we htted the individual radial prohles in Fig. ^ and obtained 
acceptable hts for all energy bands. The best £t parameters are summarized in Table ^ 
We can compare this with the results from the Einstein IPG; that is, the detection of the 
central excess emission above a /3-model (Pavid et ah 19^) . Since ASCA has a much 
poorer angular resolution than the Einstein IPG, the central excess in the IPG data may be 
unresolved in the ASCA data. Thus we obtain a single /3-model with core radius smaller 
than 1'.6, consistent with the results obtained from IPG data in the r > 1'.5 region ( Pavid 
et ah 1990|) . In fact, the GIS radial prohles are well reproduced by convolving the ROSAT 
PSPG image with the ASCA response (see § 5.3). The important result here is that all the 
radial brightness prohles in diherent energy bands show an almost identical shape, as can be 
seen from Table |l|. This suggests that the central excess emission may be present in higher 
energy bands as well as in the soft energy band. We will come back to this issue later. 


5.2. ROSAT Data and Analysis 

As described in § 4, to obtain the temperature structure in the Hydra-A cluster, we 
need to ht X-ray image and spectra simultaneously. However, the brightness prohle obtained 
from the GIS data alone would be highly model dependent because of the relatively wide 
PSFs. An alternative preferable way is to utilize the ROSAT image itself as the model of 
the brightness prohle, as has been done in Markevitch et ah (1996). We thus analyzed the 
data of the ROSAT RSRC. 

The RGS'HTPSPG observed the Hydra-A cluster from 1992 Nov 8 to 1992 Nov 9 with 
a total on-target exposure time of 18 ksec and the pointing position of 9^18”^05^, —12°06'00" 
(J2000). We obtained the processed data of the PSPG from the archival data base provided 
by the ROSAT Guest Observer Facility in NASA Goddard Space Flight Genter. For the 
data reduction we followed the procedure described by Snowden et al (1994) and used their 
software ( ^nowden 1995|) . 
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From the PSPC data, the valid time is selected so that the Master Veto count rate is 
always smaller than 170 counts s“^ to eliminate the data with anomalously high particle 
background rates (Snowden et ah 1992| ; [Plucinsky et ah 1993|) . Using Snowden’s software 
we can estimate the count rate of the four background components, the particle background 
(PB), after-pulse events (AP), scattered solar X-ray background (SB), and long-term 
enhancements (LTE). After subtracting the background, the PSPC image was corrected 
using the exposure map which represents the effective area and the exposure time as a 
function of the position in the sky coordinates. We thus obtained the background-subtracted 
flat-held images in the 0.14-2.04 keV band (Fig. |^). In the PSPC data. X-ray flux from the 
cluster is detected out to ~20 arcmin. 


The X-ray image taken by the PSPC exhibits very good circular symmetry. To 
express the X-ray brightness profile quantitatively, we performed model fitting to the 
radial brightness profiles. We assumed that the PSPC point-source responses are ideally 
point-like; this is true if we do not discuss structures smaller than ~1'. Firstly we fit the 
PSPC radial profile with a single-/? model, but the fit was quite poor. As a next step, we 
applied a sum of two /^-models to fit the profile and derived good fit as shown in Fig. 

The best fit parameters of the two /^-models are {vc, /?) = (lh36, 2.5) for the narrower 
component, and (^c/?) = (1'.80, 0.68) for the wider component. In the narrower component 
the intensity is reduced by half at r = O'.44. The normalization ratio of the two /5-models, 
/ = (normalization of the narrow /3-model) / (that of the wide /3-model), is 4.9. This 
reconfirms the existence of the central excess emission found with the Einstein IPC. 


5.3. Combined ASCA and ROSAT Data 

In the soft X-ray band where ROSAT and ASCA are both sensitive, the ASCA radial 
profile also must be fitted with the same double /3-model which fit the PSPC radial profile. 
Thus, we fit the CIS radial profile in 0.5-2.5 keV with the double /3-model, in which the 
two sets of /3 parameters and core radii (rc), and the ratio between the two normalizations 
(/) were fixed at the best fit values derived in § 5.2. We derived a good fit and the obtained 
minimum chi-square value was 20.5/(i/=19). 

In the higher energy bands, does the radial profile require a central excess? In the 
previous section, we suggested that the similar radial profiles in different energy band could 
be the evidence of the existence of the central excess emission in the higher energy bands 
as well. Using the double /3 model, we have checked if the higher energy bands also show 
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the central excess or not. In this case, we let the normalizations for the two /3-model 
components to vary independently, while the two sets of (3 parameters and core radii were 
hxed to the values which were obtained in § 5.2. The htting results are illustrated in Fig. 

In all energy ranges, the normalization ratios (/) obtained are consistent with 4.9 which 
is the best £t value obtained with the PSPC radial profile in § 5.2. Moreover, when we 
hxed the normalization of the narrow /3-component to be 0, the £t became unacceptable; 
the minimum values are 101, 133, 81, 48, and 30, with z/ = 19 for the 0.5-1.5 keV, 
1.5-2.5 keV, 2.5-4 keV, 4-6 keV, and 6-10 keV, respectively. These results imply that the 
narrower of the two /3-model components is dehnitely needed in individual radial profiles, 
not only in the soft X-ray band but also in the harder energy ranges. Therefore, we can 
conclude that the central excess emission is most likely to present in the higher energy band 
above 4 keV as well as in the soft X-ray band. 


6. TEMPERATURE AND ABUNDANCE STRUCTURE 

6.1. Temperature and Abundance Profile 

In the last section, we studied the radial brightness prohles in different energy ranges. 
In this section, we describe the temperature and abundance prohle based on the spatially 
sorted energy spectra. Using the method described in § 4, we performed simultaneous 
htting to the hve GIS annular spectra which were htted individually in § 3, fully accounting 
for the contributions of the X-ray hux from the sky to the focal plane. The image response 
matrices used for this analysis have the dimension of 26 x 5 for each of the 201 energy bins. 
Twenty-six model regions consist of 8 annuli each O'.25 wide for r = 0 — 2' and 18 annuli 
each 1' wide for r = 2' — 20'. 

As the model brightness prohle, we employed the PSPC radial brightness prohle 
obtained in § 5. We assumed that the PSPC image represents the surface brightness 
prohle in the energy range below 2 keV, and converted those surface brightness to emission 
measures according to the plasma emission code for a given temperature and abundance. 
The projected prohles of the temperature and abundance were assumed to be constant 
within the individual 5 annular regions of 0 — 1'.5, 1'.5 — 3', 3' — 5', 5' — 10', and 10' — 20', 
which correspond to the same hve annular regions used in extracting the GIS spectra in 
detector coordinates (Table 1). The hydrogen column density is assumed to be constant 
over the entire cluster. Therefore, the free parameters are 5 temperatures, 5 abundances. 
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the hydrogen column density and the overall normalization. 

The best £t parameters are summarized in Table ^ and the temperature and abundance 
prohles are illustrated in Fig. |T^by crosses. Except for the central region, there seems to be 
an indication of gradual temperature decrease towards the outside, but it is not statistically 
signihcant. In the abundance distribution, only the central value was well constrained while 
only upper limits were obtained for the outer regions. However, the central abundance value 
is larger than all other best-£t values in the outer regions. This suggests there is a central 
concentration in the abundance distribution. In order to study the difference between the 
central 1'.5 region and the outside region more clearly, we again htted the hve GIS spectra 
simultaneously assuming a common temperature and abundance outside 1'.5. The results 
are summarized in Table ^ and illustrated in Fig. 0 with diamonds. The central 1'.5 
abundance value is signihcantly larger than that of the outside, implying there is a central 
concentration in the abundance distribution. 


6.2. Central Cool Region 

Since the radiative cooling time in the central region of the Hydra-A cluster is shorter 
than the Hubble time (2 x 10^° yr), the presence of a cooling flow is expected as discussed 
by David et ah (1990). If there is no significant heat input, the thermal instability occurs 
in the central densest part of the cluster gas and a rapid radiative-cooling process would 
lead to the formation of a cooler gas phase surrounded by the hotter ICM. As presented 
in the previous sections, the ASCA spectral analysis showed that the X-ray emission from 
the central 1'.5 radius region can be attributed to the 3.4 keV gas and there is no clear 
evidence of a central cooler gas component. However, the GIS and SIS efhciencies decrease 
so rapidly below 1 keV that ASCA is not very sensitive to emission from gas cooler than 
~1 keV. The i^OFATPSPG is sensitive down to ~0.1 keV and thus these data would be 
useful to investigate whether there is an additional central cooler gas component. In the 
following, we describe the results from the PSPG spectral htting and joint analysis of the 
PSPG and GIS spectra. 

We analyzed the iiOFAT PSPG data hrst. From the same data set selected in §5.2, we 
accumulated the pulse-height spectrum in the central 1'.5 region. Each X-ray event was 
weighted using the vignetting function so that the X-ray image became flat. The original 
256 channel bins are summed up into 7 energy bands named RIL through R7 dehned in 
Snowden et ah (1994). For the model fitting, we used only the 6 bins of R2 through R7 to 
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avoid the after-pulse event contamination which could be contributing to the RIL band. 
The background spectrum was made from the annular region of r = 36' — 46' where no 
apparent point source was found. 

We htted the PSPC spectrum thus obtained with a single temperature RS-model as 
well as the MEKA-model, with the response matrix named ”pspcb_gain2_256.rsp” publicly 
available from NASA Goddard Space Flight Center. The best-£t temperature, abundance, 
and hydrogen column density and their 90% confidence errors for one parameter of interest 
obtained with the RS-model were 2.27 (2.03-2.63) keV, 0.30 (0.21-0.42) solar, and 3.79 
(3.59-3.99) xl0^° cm“^, respectively. MEKA-model gave 2.25 (2.00-2.55) keV, 0.33 
(0.23-0.45) solar, and 3.78 (3.58-3.98) cm“^. The derived temperatures are signihcantly 
lower than that obtained from the GIS data for the corresponding central region. This 
temperature discrepancy between the GIS and the PSPC spectra indicates that the emission 
consists of multiple temperatures. This is not surprising, because, even if the central 
cluster region is hlled with cool plasma, there must be an inevitable contribution from the 
foreground and background regions of the cluster where the temperature is ~3.4 keV. 

As a next step, we performed simultaneous fitting of the PSPC and GIS spectra. 
Because of the broad XRT PSF, the derived GIS spectrum within 1'.5 radius consists of 
X-ray flux not only from the corresponding sky region but also from outer regions of the 
sky. However, the contamination from outside the 0 — 1'.5 sky region does not significantly 
affect the results. This is clear from the fact that the annular spectral analysis performed 
in § 3 and the simultaneously htting of the hve GIS spectra in § 6.1 gave the same results 
for the central 1'.5 region. Therefore, we used only the central 0 — 1'.5 GIS spectrum for 
the joint analysis with the PSPC spectrum. The effective area for the GIS spectrum was 
calculated from the PSPC image within 1'.5 convolved with the GIS-I-XRT PSFs. That 
effective area represents the contribution from the central 1'.5 sky region to the 1'.5 radius 
detector region. 

We found that the single temperature model can not explain the PSPC and GIS 
spectrum simultaneously. A single temperature RS-model and MEKA-model gave minimum 
values of 89.9 and 92.4, respectively, for 60 d.o.f. There must be at least two different 
temperature components. We attempted the following two types of models to fit the 
PSPC and GIS spectra simultaneously. Model 1 is a two temperature model, that is, two 
different temperature components coexist within 1'.5. The fitting parameters are the two 
temperatures, their emission measures, common abundances, and the galactic absorption. 
We allow the emission measures of the hot component for the PSPC and GIS spectra 
to vary independently, but constrain all other parameters to be equal. We employed the 
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RS-model and MEKA-model as the plasma emission code. We derived acceptable fits 


in both cases. The best fit model is shown in Fig |^, and the derived parameters are 
snmmarized in Table The fraction of the cool component, the cool component percentage 
in the emission measnre (EM), dehned as Fcooi = E M (Tcooi) / E M{r < 1'.5), was ~2%. 
Fig. |^(a) shows the contours in the Fcooi — Tcooi space. 


Model 2 is a projected cooling flow model which is expressed as 
CF{M,Tmax, Train) + PiTmax), where CF is the cooling flow model by Mushotzky 
& Szymkowiak (1988), and P represents the additional isothermal component. CF is a 
function of the mass deposition rate M, the maximum temperature Tmax from which the 
gas cools, and the minimum temperature to which the gas cools. The temperature is 
distributed continuously from Tmax to Tmin and each temperature component has emission 
measure which is inversely proportional to the total emissivity of that temperature. The 
additional isothermal component must have the same temperature as Tmax- In this model 
2, we also applied two different emission codes, RS-model and MEKA-model (CFLOW and 
MKCFLOW XSPEC model, respectively, for CF component) and derived acceptable hts 
(Table |^). The chi-square contour map for M vs Tmax is shown in Fig[^(b). The estimated 
mass deposition rate is 60±30 Mq yr“^, while the fraction of the cooling flow component 
to the total emission measure within the 1'.5 region {Fcf = EM{CF)/EM{r < 1^5)) 
was ~8 %. With both model 1 and 2, the obtained X-ray luminosity within 1'.5 region 
was 2.6 X 10"^^ ergs s“^ in the 0.5-4.5 keV band, and is primarily the emission of the hot 
non-cooled component with the temperature of ~3.5 keV. The 0.5-4.5 keV luminosity of 
the cool component of model 1 and the cooling flow component of model 2 are 6.4 x 10"^^ 
and 2.8 x 10^^ ergs s“^, respectively. The fraction of the hot component is thus far larger 
than could be explained in terms of the projected foreground/background emission from 
the off-center regions. 


7. SUMMARY AND DISCUSSION 

The energy-dependent flux-mixing effect due to the XRT PSF makes it quite difficult to 
analyze the data of extended sources observed with ASCA. In order to perform simultaneous 
fitting of the X-ray images and spectra, we have calibrated the XRT PSF using the data 
from Cyg X-1, and have fully taken it into account in the data analysis. 

The X-ray radial brightness prohle obtained by the ROSAT PSPC exhibits a central 
excess above a single /3-model within ~1'.5, as does the radial prohle from the Einstein 
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IPC data. With ASCA^ we observed the X-ray surface brightness in higher energy bands 
up to ~10 keV for the hrst time. Using the newly developed analysis technique, we £t the 
radial brightness prohles and found that there is no clear difference among the profiles in 
different energy bands. This suggests that the central excess emission found in the soft 
X-ray images from Einstein and ROSAT also exists in the higher energy band of 4-10 keV. 
In particular, we successfully reproduced the GIS radial brightness prohle as a sum of two 
/3-models. The narrower of the two (3 components, thought to represent the central excess 
emission, is required not only below 2 keV but also in harder energy bands, up to 10 keV, 
assuming a double /3 prohle in all energy bands. 

The simultaneous fitting of the five annular spectra taken from the GIS data gave the 
radial prohles of the temperature and metal abundance (see Fig. 0 ). The obtained overall 
temperature structure is consistent with being isothermal; this result is also consistent 
with the fact that all the radial prohles are very similar. However, the PSPG spectra 
accumulated from the region within 1'.5, where the surface brightness prohle begin to 
deviate from the /3-model, gives a signihcantly lower temperature than that obtained from 
the GIS spectrum. This means that there must be an additional cool component at the 
cluster center. Therefore we jointly analyzed the GIS and PSPG data and successfully ht 
both spectra simultaneously with the two temperature model as well as the cooling how 
model (Fig |^. The cooling how model gives the mass deposition rate of 60 ± 30 Mq yr“^, 
an order of magnitude smaller than the 600 Mq yr“^ estimated from the Einstein data by 
David et ah (1990). If we assumed that all the hux coming from the central 1'.5 region is 
originated from cooling how, we would derive a value consistent with 600 Mq yr“^, using 
the formula of M = 2fimpLx/5kT and the bolometric luminosity within 1'.5 of La; ~ 4 x 10^^ 
ergs s“^. However, as we showed in § 6, the central 1'.5 region can not be entirely cooled. 
More than 90% of the total emission measure consists of the hot non-cooled component 
with the temperature of ~3.5 keV. Therefore, the central excess in the X-ray brightness 
prohle can not be formed only by the cool component. 


Since the central region representing the excess emission mostly occupied by the same 
hot IGM component which permeates the rest of the cluster, we interpret the central excess 
emission as an evidence for gravitational potential structure. The potential structure can be 
interpreted as consisting of two distinct components; a large scale cluster component and 
a central compact component attributed to the cD galaxy. In previous investigations, such 
a dual potential structure has been suggested (P?homas, Fabian, &: Nulsen 1987|; |NuIsen 
&; Bohringer 1995| ) or assumed (e.g. [Stewart et ah 1984| ). The hrst direct observational 
evidence of the additional potential dimple around NGG 1399 in the Fornax cluster was 
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found from the ASCA observation ( Ikebe et al. 1996|) . In the case of the Hydra-A cluster, 
there also must be a central potential dimple around the cD galaxy which primarily causes 
the central excess in the brightness prohle. 


Thus far, the central excess brightness seen in many clusters has been interpreted 
mainly as due to the central temperature decrease of gas, thus providing a basis for the 
cooling flow. However, our results clearly reveal that the central excess brightness is at 
least partially caused by the dual potential structure around the cD galaxy. Therefore, the 
cooling flow rate derived from the central excess brightness can be grossly over-estimated, 
as in the case of the Hydra-A cluster, if the dual potential structure is ignored. 


The metal abundance distribution is also a very important subject which is strongly 
related to the cluster evolution scenario. A detailed measurement has become possible for 
the hrst time using ASCA. In the Hydra-A cluster, we found an indication of a central 
concentration in the metal abundance distribution (Fig. 0 ). ASCA has clearly detected 
the central concentration in the heavy element abundance distributions in the nearest two 
clusters, Virgo and Centaurus (|Matsumoto et al. 1996| ; [Fukazawa et al. 19941) , which are 
the hrm conhrmation of the previous results by Ginga ([Koyama, Takano fc Tawara 1990|) 
and ROSAT ( [Alien &: Fabian 1994|) , respectively. A similar feature has also been discovered 
in the poor cluster AWM7 ( Xu et al. 1996 ). The central abundance concentration may 
be caused by a large contribution by metal-enriched ISM of the cD galaxies, which may 
have not suffered a ram-pressure-stripping process in the cluster evolution ([Tamura et al. 


19961) , because the cD galaxy is sitting in the bottom of the gravitational potential well. 


We speculate that all the clusters showing central excess emission in the brightness prohle 
have central concentration of heavy elements. 


Based on the results from the simultaneous htting of the GIS and PSPC spectra with 
the two temperature model (model 1 in § 6.2), we calculated the mass of the cool component 
gas. We assumed that the cool component coexists with the hot component within the 1'.5 
region, and the local pressure balance is achieved between the cool component and the hot 
component as UcooiTcooi = nhotThot- The mass of the cool component gas is estimated to be 
Mgas ~ 6 X 10® Mq; and its hlling factor, that is, the volume fraction of the cool component, 
is ~6 xl0“^. On the other hand, if the cool component is concentrated at the center, and 
the pressure equilibrium is achieved at the boundary of the cool component and the hot 
component, the cool component would be distributed out to ~5 kpc and would have a 
total mass of Mgas ~ 2 x 10® M©. These values are about 0.1-0.4% of the stellar mass of 
the cD galaxy 3C218 = Hydra-A, estimated to be 1.5 x 10^‘^Mq based on the assumption 
of M/Lb = Q{M/Lb)q and logL^ = 11.4. This is comparable to that in other non-cD 
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ellipticals ([bbrman, Jones fc Tucker 1985|) . 
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snpporting to bnild the data analysis tools. We also thank all members of the ASCA 
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for important comments and helping PSF calibration, and to Joel Bregman for referee’s 
comments. Y.I. acknowledges snpport from the Special Researchers’ Basic Science Program 
and thanks Makoto Hattori for helpfnl discnssion. M.H. acknowledges snpport from the 
Research Fellowships of the Japan Society for the Promotion of Science for Yonng Scientists. 


A. PSF CALIBRATION 


For the analysis of the extended sonrces, adeqnate knowledge of the XRT PSF is 
reqnired. However, neither the calibration experiment performed on the gronnd nor the 
ray-tracing code conld reprodnce the PSFs with snfficient accnracy, and in-orbit calibration 
was reqnired. As we described in § 4, the PSF has largely extended ontskirts and depends 
on the X-ray energies. In addition to these characteristics, we shonld note the difhcnlties 
of the PSF modeling cansed by its lack of cylindrical symmetry and strong dependency on 
the position on the focal plane. Accordingly we have developed the method to reprodnce 
the PSFs by interpolating among a set of data from Cyg X-1 taken at varions positions on 
the focal plane (fTakahashi et ah 1995| ; |Ikebe 19^) . Cyg X-1 was selected becanse it is one 
of the brightest point X-ray sonrces with a hard spectrnm. With this method we can only 
reprodnce the XRT-I-GIS PSFs. Becanse any point source bright enough to yield a sufficient 
number of signal photons within realistic observing times would cause considerable event 
pile-up in the SIS data, the XRTj-SIS PSFs is difficult to obtain. 

The observations of Cyg X-1 used to produce the XRT-I-GIS PSFs were performed 
in November 1993, November 1994 and May 1995 at several offset positions. Due to the 
5' relative misalignment between the two XRTs equipped the two CIS sensors, GIS-S2 


and GIS-S3, we obtained 14 sample positions on the focal plane (Fig |T^. Assuming that 
the eight quadrants of the two XRTs are identical and perfectly symmetric, we can have 
practically 92 PSF data sets. We have divided each pointing data set into 8 energy bands, 
0.5-1 keV, 1-2 keV, 2-3 keV, 3-4 keV, 4-5 keV, 5-6 keV, 6-8 keV and 8-12 keV, yielding a 
PSF image data base composed of 8 x 92 PSF images. When we produce a PSF at a given 
position and energy, we select the nearest 2-8 sample positions in Fig. O and the nearest 2 
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energy bands, and linearly interpolate them in the position and energy space. 

We examined the accuracy of the generated PSFs by comparing the synthesized PSF 
with the observed data of 3C273, in terms of their radial profiles centered on the peak. At 
the energies above 2 keV, the radial prohles agree reasonably. On the other hand, below 
2 keV band, the PSFs generated from Cyg X-1 data have larger HPD by ~ 0.'2 and are 
signihcantly wider than those of 3C273. This may be caused by the intrinsic extension in 
the X-ray emission from Cyg X-1 due to the X-ray scattering by interstellar dust grains 
( [Mitsuda et ah 1990|) . Therefore we improved the PSFs below 2 keV as follows. For the 
1-2 keV band PSF images in the data base, we combined the outskirts region (r > 6') of 
the 1-2 keV band with the core (r < 6') of the 2-3 keV band images convolved with a 
Gaussian of cr = 0M25. As the 0.5-1 keV band images, we used the 2-3 keV band images 
which is smoothed with a Gaussian of cr = 07625. 


Fig [Tj shows the comparison of the PSFs thus obtained with those of the data from 
3G273, in terms of their radial profiles centered on the peak. Based on studies using 
Monte-Garlo simulation, we introduced a systematic error 5% for the model predictions in 
the analysis in this paper. The systematic error of the produced PSFs may be caused by 
the inaccuracy of the assumption that the eight quadrants of the two XRTs are identical 
and symmetric. Futher, additional observation of Gyg X-1 would reduce the systematic 
error in the interpolation procedure on the focal plane. 
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Fig. 1.— The X-ray intensity contour map of the Hydra-A cluster obtained with the GIS. 
The background subtracted image was convolved with a Gaussian of 1' FWHM. The contours 
are logarithmic where each step corresponds to a multiplicative factor of 2. 

Fig. 2.— The GIS (a) and SIS (b) spectra extracted from the central 1.5' radius region 
htted with a thermal bremsstrahlung plus three emission lines modihed by the photo-electric 
absorption. Grosses are the data and the histograms are the model. For the display, SIS-Sl 
data and its model are multiplied by 0.1. 

Fig. 3.— The X-ray flux contributions from the sky annular regions to the individual 
annular spectra of the GIS are illustrated. We simulated the GIS data of an isothermal 
symmetric cluster having the radial brightness prohle shown in Fig. 8. The temperature, 
heavy element abundance, and hydrogen column density are assumed to be 4 keV, 0.3 solar, 
and 5 x 10^'^ cm“^, respectively. On the horizontal axis the radius ranges of extracted spectra 
on the detector are plotted. The vertical axis shows the contribution fraction from different 
sky regions, illustrated with different patterns. 

Fig. 4.— The results of the annular spectral fits with a thermal bremsstrahlung plus three 
emission lines. Temperature (a) and equivalent width of Fe-K lines (b) are illustrated. 

Fig. 5.— The results of the annular spectral hts with RS-model. Temperature (a) and 
abundance (b) are illustrated. 

Fig. 6.— The background-subtracted radial brightness prohles in different energy ranges 
obtained by the GIS. 

Fig. 7.— The X-ray intensity contour map obtained with the ROSAT PSPG. The image 
was smoothed with a Gaussian of 1' FWHM. The contours are logarithmic where each step 
corresponds to a multiplicative factor of 2. 
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Fig. 8.— The background-subtracted PSPC radial profile fitted with the double /5-model. 
Open circles are the data and histograms are the models. 

Fig. 9.— The best £t normalization ratios between narrow f3 and wide (3 components. The 
errors represent 90 % confidence limit. 

Fig. 10.— The radial profiles of the temperature (a) and abundance (b) derived from the 
simultaneous htting of the five GIS annular spectra. 

Fig. 11.— Simultaneous fitting of the PSPC and GIS spectra extracted from the central 1.5' 
radius region. Filled circles are the PSPC data and crosses are the GIS data. Histograms 
are the best fit two temperature model. 

Fig. 12.— The contour map (a) for the temperature of the cool component (Tcooi) vs the 
fraction of the cool component {Fcooi, see text) and (b) for the mass deposition rate (M) vs 
maximum temperature (Tmax)- 

Fig. 13.— The positions on the focal plane where the Cyg X-1 is observed (hlled circles) 
and the positions which are identical to those assuming symmetry (open circles). 

Fig. 14.— Comparison of the radial brightness profiles of the 3C273 data with those of the 
PSFs that we used in our analysis. The PSFs are reproduced by interpolating among the 
actual observed data of Cyg X-1 or their modihed images (see text). The vertical axis shows 
the radial prohle ratios (data/PSF) averaged over eight data sets of 3C237. 



Table 1. The results of the GIS radial prohle htting with a single /3-model. 


energy range(keV) 

rc( arcmin) 

/9 


0.5-1.5 

0.97 (0.71-1.25) 

0.62 (0.58-0.66) 

10.5/17 

1.5-2.5 

0.96 (0.70-1.30) 

0.63 (0.59-0.68) 

15.1/17 

2.5-4 

1.14 (0.69-1.35) 

0.67 (0.58-0.68) 

7.0/17 

4-6 

1.54 (1.10-2.25) 

0.72 (0.64-0.87) 

5.6/17 

6-10 

1.20 (0.58-2.05) 

0.67 (0.56-0.87) 

7.3/17 


Note. — The errors represent 90% conhdence limit for two 
parameters of interest. 


Table 2. Simultaneous htting of the hve GIS annular spectra 


Radius(arcmin) 

0-1.5 

1.5-3 

3-5 

5-10 

10-20 

Temp er at ure (ke V) 

o HI +0.20 

^•+*0-0.81 

4.0811:?^ 

q co+O.OS 

c\ ^^+0.75 
' '-0.60 

Abundance (solar) 

n c;7+0.25 

U-O'-0.21 

0.17 (<0.55) 

0.015 (<0.74) 

0.25 (<0.71) 

0.25 (<0.94) 

Ah (x 10^° cm“^) 



1.7 (< 3.9) 






237.1/283 




Note. 


The errors represent 90% conhdence limit for one parameter of interest. 











Table 3. Obtained temperatures and abundances in the central 1'.5 region and the outside 


Radius(arcmin) 

0-1.5 

1.5-20 

Temperature(keV) 

3.55±0.18 

3.9il[;:i 

Abundance (solar) 

0.56l°:is 

0.14±0.11 

Nil (xl0^° cm“^) 

1.2 ( 

< 3.2) 


247.3/289 


Note. — The errors represent 90% 


conhdence limit for one parameter of interest. 







Table 4. Simultaneous fitting of the PSPC and GIS spectra with a two temperature 

model 


plasma 

T 

cool 

That 

EM (Tcooz) 

EM{Thot) 

abundance 

Nh 


code 

(keV) 

(keV) 

(cm“^) 

(cm“^) 

(solar) 

(cm“^) 


RS 

0.67 

3.45 

3.9 X lO®® 

2.4 X 10^^^ 

0.51 

3.4 X lO^o 

74.15/58 

MEKA 

0.50 

3.41 

4.1 X lO®® 

2.4 X 10^^^ 

0.56 

3.4 X lO^o 

73.97/58 


Note. — Both the PSPC and GIS spectra were extracted from the central 1'.5 region. 
The htting model is expressed with sum of two RS-model or MEKA-model modihed by 
photo-electric absorption. 


Table 5. Simultaneous fitting of the PSPC and GIS spectra with a cooling flow model 


plasma 

M 

That 

EM{CF) 

EM{Thot) 

abundance 

iVu 


code 

(Mo yr-i) 

(keV) 

(cm“^) 

(cm“^) 

(solar) 

(cm“^) 


RS 

57 

3.57 

2.0 X 10^^^ 

2.3 X 10*^^ 

0.50 

3.3 xlO^o 

73.85/59 

MEKA 

62 

3.56 

2.2 X 10^^^ 

2.2 X 10*^^ 

0.51 

3.4 xlO^o 

73.72/59 
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